ac-driven atomic quantum motor 
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We invent an ac-driven quantum motor consisting of two different, interacting ultracold atoms 
placed into a ring-shaped optical lattice and submerged in a pulsating magnetic field. While the 
first atom carries a current, the second one serves as a quantum starter. For fixed zero-momentum 
initial conditions the asymptotic carrier velocity converges to a unique non-zero value. We also 
demonstrate that this quantum motor performs work against a constant load. 
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Linear or rotational motion presents the basic work- 
ing principle powering all sorts of machines. For nearly 
two centuries, since the invention of the first electrical 
motor [H , the ever continuing miniaturization of devices 
has profound consequences for several branches of sci- 
ence, industry, and everyday life. This process has al- 
ready passed the scale of micrometers |2| and has en- 
tered the realm of the world of nanoscale [3( . Bioinspired 
devices such as chemical or light driven synthetic molec- 
ular motors identify just one of those recent successes 
[3| . While the operational description of such molecular 
motors mainly rests on classical concepts, much less is 
known for operational schemes that are fully quantum 
mechanical in nature. An ideal resource for the latter 
possibility is the dynamics of cold atoms that are posi- 
tioned in optical potentials [5j]. 

With this work, we put forward a setup for a quan- 
tum motor which consists of two species of interacting, 
distinguishable quantum particles that are loaded into 
a ring-shaped optical potential. The blueprint for such 
an underlying ring-shaped one-dimensional optical lattice 
has been proposed recently [6| and a first experimental 
realization has been reported in Q. Here, we employ 
this setup to devise an engine which works as a genuine 
ac-driven quantum motor. 

ac-quantum motor. — Figure 1 outlines our device. The 
ring-shaped optical potential, which results either from 
the interference of a Laguerre-Gauss (LG) laser beam 
with a plane wave Q or, alternatively, of two collinear LG 
beams with different frequencies 0] is capable of trapping 
two interacting atoms. One of the atoms, termed "car- 
rier" , c, is driven by an external field, while the other 
atom, termed "starter" , s, interacts locally via elastic s- 
wave collisions with the "carrier" , but remains unaffected 
by the driving field [§] . Two possible setups that come to 
mind are (i) A neutral "starter" and an ionized "carrier" , 
a suitable driving field can be implemented in a way typ- 
ically done for electrons placed in a conducting ring, i.e., 
by a time-dependent magnetic flux threading the lattice 
[9(. (ii) A spinless "starter" and a "carrier" -atom with a 
non-zero spin which is driven by a time-dependent cone- 
shaped magnetic field of an Ioffe-Pritchard trap @, [l(| • 

We next assume that both atoms are loaded into the 
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FIG. 1: (color online) Atomic quantum motor: Two different 
ultracold atoms are loaded into a ring-shaped optical lattice. 
Both atoms interact locally with each other, while only one 
carrier (the one with an arrow) is magnetically powered. 



lowest energyband of a deep, ring-shaped optical po- 
tential with L lattice sites and the lattice constant d. 
The time-dependent homogeneous vector potential A(t) 
does not induce any appreciable transitions between the 
ground band and the excited band(s). 

The ac-driven, total Hamiltonian fltot of the motor 



H tot = H c (t) + H s + H i 



int 



(i) 



is composed of the time-dependent Hamiltonian H c (t) for 
the "carrier" 



H c (t) = ^^ e ^W|? c + l)(/ c | +H.c. 

and for the "starter" H s , respectively, i.e., 
l \ 



H.c. 



(2) 



(3) 



Here, J c and J s are the corresponding hopping strengths 
which are functions of the atom masses and the optical 
potential depth 11]. The salient carrier-starter (on-site)- 
interaction reads 

L 

H mt = W Y, Kjl.\Ic)(Ic\®\L)(Ib\, (4) 
i c ,i s =i 

where W denotes the interaction strength. Throughout 
the remaining we use periodic boundary conditions; i.e., 
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\L + 1) = |1). The full system Hilbert space is spanned 
by the direct products of single particle Wannier states 
\lc) ® \h), with the dimension being TV = L 2 . The scale 
of the motor current will be measured in units of the 
maximal group velocity vq — J c d/h. 

The driving of the ac-atomic quantum motor is 
switched on at the time instant to, so that the vector 
potential assumes the form 



A(t) = x (t - t )A(t) . 



(5) 



where x(t ~ *o) is the step function, and A(t) is defined 
on the entire time axis, t G (— oo, +oo). 

dc- quantum current. — The mean "carrier" current is 
given as the speed of the motor by using the velocity op- 
erator: v c (t) = i/H[H tot (t),x e ]. With x c = J2i l c\l c ) (l c \, 
one finds v c (t) = -i(v /2) (££ =1 e lA ^ \l c + 1) (l c \ 
— H.c.) ® l s - In the quasimomcntum representation 
with \k) = ^ , exp(iKn)\n), its quantum expectation 
v c (t;t ) = (ip(t)\v c {t)\i/j(t)) reads 

L 

v c (t; t ) = vo ^2 pm c (t; to) sin (ki c + A(t)J , (6) 

ic = l 

wherein m c = 2~kI c /L is the single particle quasimo- 
mentum and where we indicated its parametric depen- 
dence on the start time to- Further, p Kl (t;to) = 
\(i>(t)\ni e ) ® |k/ c }| 2 , with K ie = 2nl s /L, is the quasi- 
momentum distribution for the carrier. 

The steady state regime of the motor can be charac- 
terized by the dc-component of the averaged velocity 



v c (t ) := lim t _ 



v c {t';t )dt'. 



(7) 



In the absence of the interaction between the particles, 
i.e., W = 0, and an initial preparation with localized 
carriers that start out with zero velocity the motor set- 
up cannot even support a transient directed current. In 
fact, this result holds for any shape of the vector potential 
A(t) [13]. This situation thus mimics a single-phase ac- 
motor: a periodically pulsating magnetic field would fail 
to put a rotor from rest into rotation, unless one applies 
an initial "push" via a starter mechanism [l3j |. In our 
setup, the role of the quantum starter is taken over by 
the non- vanishing interaction W with the second particle. 

Nonetheless, even with nonvanishing interaction, 
\W\ > 0, there exists no evident procedure for setting 
the motor into rotation. A seemingly obvious solution 
- application of a constant carrier-bias - cannot resolve 
the task. This is true because the corresponding vec- 
tor potential, As(t) = cost, induces Bloch oscillations 
only [14| . In distinct contrast, we shall employ an un- 
biased time-dependent vector potential possessing a zero 
dc-component, A(t + T) = A(t), i.e., J Q T A(r)dr = 0. 

For zero-momentum initial conditions, the unbiased 
monochromatic ac- force, A(t) = Asin(u>t), would launch 



- with equal probabilities - the system either into a clock- 
wise (rightward) or a counterclockwise rotation (leftward 
motion) [l5j]. Thus, the modus operandi as a motor re- 
quires a symmetry-breaking driving field, realized here 
with the harmonic mixing signal, 



A(t) = Ai sin(urf) + A 2 sin(2wt + 6) 



(8) 



where denotes the crucial symmetry-breaking phase 
shift. The input ([5]) knowingly may induce a non- 
vanishing nonlinear response, the so-called ratchet effect 

II EE m 

Quantum current in terms of Floquet states. — The dy- 
namics at times t > to of the time-periodic Hamiltonian 
(fT]) can be analyzed by using the Floquet formalism ■ 
The solution of the eigenproblem: U(t, to)\(f> n (t; to; k)) — 
exp(— |e n t) \4> n (t;to;k)), with the propagator U(t, t ) = 

Texp (— ^ J t * if tot (r)^ dr (T denotes the time ordering), 

provides the set of Floquet states, being time-periodic, 
i.e., with T — 2ir/ui being the driving period, \4> n {t + 
T;t ; fc)) = \(f> n (t;t ; k)). Here, k = J2 l>m {<pi \ki) b <8> |« m } c 
is the total quasimomentum of the Floquet state. Be- 
cause of the discrete translation invariance of the sys- 
tem, the total quasimomcntum is conserved during the 
time evolution, thus serving as a quantum number. Since 
-fftot is a function of the time difference t — to only, 
the quasienergies e„ are independent of to, and the Flo- 
quet states for different start times to are related by 

|<Mt;t ;fc)} = \Mt-t o ;0;k)). 

Using this relation we next decompose ip(to) in 
the complete basis of Floquet states: \ip(to)) = 

E^Li c ral^n(°;*0; fe )) = Y^lCnltniT - to;0;k)). In 

the absence of the driving, A(t) = 0, the motor set- 
up |l][3]) possesses the continuous translational symme- 
try in time. Thus, the expansion coefficients of the ini- 
tial wave- function ip(to) in the system eigenbasis know- 
ingly do not depend on time. On the contrary, eigen- 
states of a periodically driven system - the Floquet 
states - evolve in time, being locked by the external 
ac-field. The expansion of an initial wave-function over 
the Floquet eigenbasis depends on the start time to HJ, 
which determines the phase of the driving ac-field [la ]: 
c-n = c n (to) = (4>n{T — to, 0; k)\ip(to)). Substitution of 
the above decomposition into yields the result 



v c (t ) 



n=l 



V n |c„(t )| 2 , V 



1 

T 



v n {r)dT. (9) 



Here, v n denotes the quantum average velocity of the 
n-th Floquet state ©■ Because the Floquet states are 
periodic functions of the time difference t = t — to only, 
the velocities v n do not depend on to, and the depen- 
dence of the generated dc-current on the to solely stems 
from the coefficients c n (to). Since the system evolution 
is fully quantum coherent, i.e. there is no memory eras- 
ing induced by an environment, the asymptotic current 
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FIG. 3: Dispersion of the motor velocity versus the num- 
ber of lattice sites L. Here = 7r/2, and the other parameters 
are the same as in Fig. [2] 



FIG. 2: (color online) Averaged motor velocity in (|10p (in 
units of the maximal group velocity vo = J c d/h) as a func- 
tion of the phase shift in for L = 16. The (io)-averaged 
velocity (0 obtained by the direct time propagation of the ini- 
tial state up to 200T (dashed line) is compared to the asymp- 
totic dependence given by the Floquet approach @ (red solid 
line). The parameters are hu = 0.1 x J c , A\ — 0.5, A2 = 0.25, 
W = 0.2J C , J s = J c = J. 



maintains the memory of the initial condition as encoded 
in the coefficients c n (to) flBj . 

Input/Output characteristics. — The question is now, 
how can we control the motor? To answer this ques- 
tion, we used the symmetry analysis [15| which allows us 
to predict an appearance of a certain dc-current. Com- 
bining time-reversal operation and the complex conjuga- 
tion applied to (O with A(t) in the form of (JSJ), one can 
prove the (anti-) symmetric dependence of v n on for 
the Floquet states with k = 0: v n (ir — 0) = v n (Q), 
v n (— 0) = —v n (Q). Thus, the Floquet states with 
k = possess zero mean velocities at O = 0, ir. Fur- 
thermore, using a similar reasoning, one finds that the 
set of Floquet states with nonzero k can be ordered by 
the parity relation, which links eigenstates with opposite 
quasimomenta, (j> n (t;to; —k;Q) = 4> m {T — t; to; k; — 0), 
yielding v n = —v m . This implies that for a symmetric 
(in k) initial state and = 0, 7r, the contributions to 
the dc-current of Floquet states with opposite quasimo- 
menta eliminate each other. The same holds true for a 
monochromatic driving ([8]) , with A2 = [IB} . Shifting © 
away from 0, ±7r causes the decisive symmetry breaking 
and leads to the desymmetrization of the Floquet states 
with k = and consequently will violate the parity be- 
tween states with opposite signs of k. 

The motor speed depends on the initial conditions, 
which define the contributions of different Floquet states 
to the carrier velocity We restrict our analysis to the 
initial state ip(t a ) = L~ 1//2 |^ c ) (g) Y^i \ls), l c = 1,...,L, hi 
the form of the localized carrier (at l c ) and the uniformly 
"smeared" , delocalized starter. Both particles have zero 
velocities at t = to ■ The asymptotic velocity typically ex- 
hibits a strong dependence on to [lj| ■ We first discuss the 
results obtained after averaging over to, thus assigning a 



unique motor velocity value, 



= (v c (t )) ta = 1/T 



T+t 



v c (t )dt , (10) 



for fixed system parameters. 

Figure [5] depicts the dependence of the average motor 
velocity on 0. The results obtained by direct time prop- 
agation of the initial state and averaged over to (dashed 
line) are superimposed by those calculated via the Flo- 
quet formalism © (solid line). The agreement between 
the two curves is satisfactory but not perfect: This is true 
because of the sharp peaks on the asymptotic current ([9]) . 

These peaks can be associated with avoided crossings 
between two quasienergy levels 15( . These avoided cross- 
ings cause a strong current enhancement if one of the 
interacting, and transporting eigenstate overlaps signifi- 
cantly with an initial, nontransporting state of the mo- 
tor. Note also that a very narrow avoided crossing re- 
quires a very large evolution time to become resolved, 
i.e., t b s ~ h/\e a — ep\ [15j |. Our chosen evolution time 
t = 200T is not large enough to clearly resolve the dis- 
tinct resonances depicted in Fig. [21 

We further detect that the dependence of the motor 
velocity v c (to) in ([7]) on to increasingly disappears upon 
increasing the size L (not shown) . To provide a quantita- 
tive estimate, we evaluated the dispersion of the current 
© with respect to to, i.e., 



^(v c {to? (v c (t ))l) t • (11) 



Co 

As shown in Fig. [3] this dispersion decays with increasing 
size L, being rather faint for L > 16. For sizes L > 16 the 
carrier assumes an asymptotic velocity that essentially 
is independent of the initial start time to- This effect is 
caused by the presence of the starter: The carrier velocity 
is obtained as the trace over the part of the total system 
Hilbert space, associated with the starter. The starter 
dynamics mimics a dissipative, finite heat bath for the 
carrier dynamics whose effectiveness increases with both, 
the (i) the dimension of the starter subspace, i. e. the 
size L, and (ii) the strength of the interaction W. 

Load characteristics. — The analysis based on Eqs. |T]- 
131 has been for a free rotator. In order to qualify for a 
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FIG. 4: The average motor velocity for "resonance" driv- 
ing and an initial condition "localized carrier and delocal- 
ized starter" (r = 10). The parameters are W = 0.2 J c , 
J s = J c = J, Ficu = O.Uc, Ai = 0.5, A 2 = 0.25, 9 = tt/2, and 
L = 4. 



genuine motor device, the engine must be able to operate 
under an applied load. The load is introduced as the bias 
uist, being added to the vector potential A(t). All the 
information about transport properties can be extracted 
by using again the Floquet formalism, provided that the 
ac-driving and the Bloch frequencies are mutually in res- 
onance [l9(, i.e. qui — tujb, where r and q are co-prime 
integers. Figure 4 shows the dependence of the asymp- 
totic motor speed for different bias values. There are two 
remarkable features. First, the spectrum of velocities is 
symmetric around lob = 0. This follows because of the 
specific choice of the phase shift at O = tt/2. Second, 
while some regimes provide a transport velocity along 
the bias, others correspond to the up-hill motion, against 
the bias. Therefore, a stationary transport in either di- 
rection is feasible. The load characteristics exhibits a 
discontinuous, fractal structure and, in distinct contrast 
to the classical case , it cannot be approximated by a 
smooth curve. This is a direct consequence of the above 
mentioned resonance condition. 

Experimental realizations. — For an experimental real- 
ization of this quantum atom motor the following features 
should be respected: (i) In the case of the setup "carrier 
with a spin/spinless starter", the carrier should assume a 
magnetic number mp > 2 [6] , to efficiently induce the ac- 
field amplitudes A\ t i- (ii) Because in the tight-binding 
approximation the maximal amplitude of the tunneling 
is limited from above, J c < Jmax — 0.13-Eq, for 6 Li atom, 
the lattice spacing d ~ [3], and hu> = 0.1J C (used 

in the calculations), the driving frequency u> should be 
less than 2Hz. Then, the time required to launch the 
motor (i.e., to approach the asymptotic velocity value) is 
around a minute. Further focusing of the laser beam can 
decrease the lattice constant d, thereby decreasing the 
launch time to experimentally accessible coherence times 
around 10 seconds 
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Conclusions. — We studied a quantum ac-motor made 
of the two species of ultracold interacting atoms, i.e. a 
"carrier" and a "starter" , moving in a ring-shaped trap- 
ping potential. For zero-momentum initial conditions 
the asymptotic carrier velocity loses its dependence on 



the switch-on time to of the ac-drive upon increasing the 
lattice size L. A natural question that arises is, What 
about the averaged starter velocity v s l We find that the 
latter sensitively depends on the system parameters: It 
can either be very small compared to the carrier velocity 
or also larger than v c . In short, the starter can move co- 
directionally or contradirectionally to the carrier motion. 
Finally, an extension of our motor setup to several inter- 
acting bosons (i.e. a finite bosonic "heat bath") presents 
an intriguing perspective. 

We acknowledge I. V. Ponomarev for providing us with 
the illustration shown in Fig[TJ This work was supported 
by the DFG through grant HA1517/31-1 and by the Ger- 
man Excellence Initiative "Nanosystems Initiative Mu- 
nich (NIM)". 
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